Irregular Dynamics in a Solvable One-Dimensional Quantum Graph 
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We show that the quantum single particle motion on a one-dimensional line with Fülop-Tsutsui 
point interactions exhibits characteristics usually associated with nonintegrable systems both in 
bound state level statistics and scattering amplitudes. We argue that this is a reflection of the 
underlying stochastic dynamics which persists in classical domain. 
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The advancement in nano-engineering in the last 
decade has brought novel incentives to the study of 
low-dimcnsional quantum systems with geometrically de- 
signed forms that have no counterpart in nature. The 
quantum graph, which is a genèric one-dimensional 



model of nano-device composed of quantum wires, repre- 
sents one of such systems 0, Q . The interest to the quan- 
tum graph is enhanced with its possibility to emulate the 
two-dimensional system of quantum billiard |3|], whose 
solution has required rather extensive numerical treat- 
ments. It is therefore quite appropriate, at this point, to 
investigate genèric aspects and general features of quan- 
tum graphs ahead of detailed studies of specific models 
of nano-devices. 

In a parallel development, quantum graphs have been 
used as a tractable model for the study of quantum chaos, 
or the irregular aspects of quantum dynamics occurring 
as quantum manifestations of classically chaotic systems 
0, 0. Naturally, it is expected that random quantum 
graphs, which are complex networks of quantum lincs, 
would result in the universal quantum fluctuation that 
has been associated to the quantum chaotic dynamics 
0,0- It has been revealed, however, that no real complex 
network is required for the irregular quantum dynamics 
to present itself. A very simple version of the quantum 
graph, the star graph, which is a quantum graph with 
many lincs conncctcd at a single node, has been shown to 
display the characteristics associated to partial quantum 
chaos in an analytical semiclassical study A natural 
question to be asked is whether we can further simplify 
the model of quantum chaos to the point of solvability. 

In this article, we consider one of the simplest possible 
quantum graph which is made up solely of nodes with two 
connected lines with the special property for the nodes 
called scale invariance. The resulting system amounts 
to a single one-dimensional line with number of scale in- 
variant point interactions. We show that the system has 
elementary analytical scattering matrices and also ele- 
mentary analytical eigenvalue equation, yet displays full 
characteristics of irregular quantum dynamics, both in 
scattering amplitudes and in bound state level statistics. 

We discuss the implication of the results, and look into 
the apparent contradiction of the appearance of the quan- 
tum chaos in a seemingly integrable, solvable conserva- 
tive one-dimensional system. 



FIG. 1: A schematic depiction of our model system. A free 
quantum particle moves along the line on which point interac- 
tions described by (3), (4) and (6) are located at si, • ■ ■ , sjv- 
This example shows the case of A = 5. 



We consider a quantum particle, constrained to move 
on a one-dimensional line with N point-like defects 
whose locations are given by x — Sj with i = 1,2,. .N 
(FIG. 1). The Hamiltonian of the system is given, in 
appropriately rescaled unit, by 



2olx 2 



on (1) 
x £ (-oo, si) U (si, s 2 ) U ...(s/v-i, s N ) U (sjv, oo). 



The dynamics of the system is described by the 
Schròdinger equation 



Htp(x) = Ei/j(x), 



(2) 



supplemented by the U(2) connection conditions at the 
defects 0, which is conveniently specified 0] by 



(Ui - /)$(*) + iL (U % + /)$'(*) = 0, 



(3) 



where i runs as i = 1, 2, N, and Ui is a unitary matrix 
belonging to U(2) group. The boundary vectors <fr(x) 
and $'(x) are given by 



-Tp'-(Si) 



(4) 



where ip±(si) and tp'±(si) denote the limit value of ipi x ) 
and its derivative from the upper and lower regions of the 
defects Sj, x — > Sj ± 0. The constant Lq is a length scale 
introduced to account for the scale anomaly |l2j | inherent 
in one-dimensional point interaction. For technical sim- 
plicity, we assume all Ui to be identical, U = Ui. Among 
all possible U , we consider scale invariant subfamily, dis- 
covered by Fülòp and Tsutsui whose U has property 



det[U±I] = 0. 



(5) 
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This condition guarantees the equation © without any 
involvement of the scale parameter L . The Standard 
parametrization for this class of U is 



U 



cos 9 é 1 ^ sin 9 
e^sinfl -cosé» 



(6) 



This gives the connection condition which reads 



gí0 

a 

é*a.il>'_(s i )=il>' + {s i ), 
where the "strength" a is defined by 

9 

a = — cot — . 



(7) 



(8) 



The Fülòp-Tsutsui point interaction © is a less known 
subclass of one-dimensional point interaction compared 
to the Standard S potential and 5' (or e) potential, but its 
property of scale invariance comes in handy in our follow- 
ing treatment. We stress that this seemingly exòtic in- 
teraction is nevertheless realizable as a short-range limit 
of certain local potential . We first consider the scat- 
tering by a single defect. Let us, for a moment, supposc 
that there is only a single defect located at x = s. Con- 
sidering the incoming wave from x < s\ side, we assume 
the wave function to be in the forní 



ip(x) 



Ri{si)e 



-ikx 



(x < Si), 



= T x {si)e ikx (x> Si ). (9) 

We obtain the transmission and reflection amplitudes as 

2a tj. „ , > 1 - a 2 



TM) 



o Í! t> 



l + a 2 



Ri(sí) = 



2 i k s 



l + a 2 



(10) 



J 



For the scattering from x > sn side, we write 



-ikx 



(x > sí). 



= e~ ikx - R^yi™ (x<8i), (11) 

and obtain 

T{{ Si ) = T*( Si ), R'M) = -RM). (12) 

The absence of the scale parameter Lo results in the en- 
ergy independence of scattering amplitudc. 

With elementary àlgebra, we can write the scattering 
amplitudes for TV defects in the recursive forms; For the 
lcft-right amplitudes TV and Rn, we have 



Ti(si)Tjy-i(s 2 , sn) 
1 + Rl(si)R N - 1 (s 2 , ...,sn)' 

Ri(si) + R N - 1 (s 2 ,...,s N ) 
1 + Rt(si)R N -i(s 2 , ...,s N )' 

The right-left amplitude T' R' are obtained from 



T/v(si, sn) = 
Rn{si, sn) — 



(13) 



T' N (si,...,s N ) = T N (s N ,—,si), 
R N (si, ...,sn) = —R* n {sn, — , si). 



(14) 



Note the reversed ordering of Sj in right and left hand 
sides of the equations. With repeated iteration, we obtain 
explicit expressions for scattering matricesin the form 



T N (k) 
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JY 



R N (k) 



B N {k) 



D N (kY ' D N (k) 

where B-N(k) and DN(k) are defined by 



(15) 



N 



N 



N 



B N (k) =/3j2 e2ÍkSt +f3 3 J2 e2Ík(Sl· ~ Sj+Sm) +f 35 J2 e2Ík{S '~ S3+Sm ~ S " +Sp) + '''> 



i>j>m 



i>j>m>n>p 



N 



N 



N 



DN{k) = 1 + 1 ^ e 2ifc(si-Sj) + ^4 e 2ifc(s 4 -Sj+s m -s ra ) _|_ p6 e 2ik(s i -s j +s m -s n +s p -s q ) _|_ 

2> i 7>m>n i>j>m>n>p>q 



r 



(16) 
(17) 



and the abreviations 







1 



7 



2a 



J<t> 



(18) 



l + a 2 ' ' l + a' 
are used. The sum runs over all indices in the range 
between 1 and N with the specified constraint, and the 
numerator contains terms up to the order of [3^ N ^ 2 ^ where 
the exponent signifies the integer part of N/2. For given 
N, there are nCi terms with order (3 l , and the scattering 
matrices are the multi-periodic oscillatory functions with 
2-Aí— i freqüències. 



Along with scattering, we can also consider the bound 
spectra by limiting the system to finite line of size L > 
sn- One of the easiest way is to impose Dirichlet bound- 
ary conditions at x = —L and x — L, ip(L) — tjj(—L) = 0. 
This leads, for the case of <j> = 0, to the eigenvalue equa- 
tion 

(R N (k) - e - 2ïfcL ) (R' N (k) -e~ 2lkL ) = T N (k)T' N (k).(19) 
Explicito calculation again yields the form 



3 



JY 



N 



N 



sin2fcL + /3 ^ sin 2fcs; + f3 2 J^sin 2k(sj-Si +L) + /3 3 y^sin 2k(s m - Sj + Sj) 

i i>j i>j>m 

N N 

+j3 i ^2 sm2k(s n -s m +Sj-Si+L) + /3 5 V] sin2£;(é 



2>j >m>n 



i>j >m>n>p 

1 



•• = 0, 



(20) 
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FIG. 2: Transmission probability as the function of inci- 
dent momentum. A common strengths parameter 01 = 3/2 is 
adopted. The location of the points are chosen to be si = 1, 
S2 = si + \/2, S3 = S2 + \/3, s 4 = S3 + V5 and s 5 = s 4 + \/7, 
se = «5 + and S7 = se + a/Ï3- 



which is a bound state counterpart of (|15|) . 

In order to reveal the physical content of the scatter- 
ing matrices and the spectral function 1)19(1 . we plot 
|T/v| 2 as the function of incident momentum k for value 
of a = 3/2 in FIG. 2. The number of point defects is 
set to be N = 3, TV = 5 and N = 7. The angle tf> is 
set to be zero for all cases. The locations Sj are set to 
be the sum of square root of primes; s, = J2]=i \/Pj-> 
where pj stands for the j-th smallest prime with con- 
vention p\ = 1, p 2 = 2. These vàlues are selected to 
guarantee the incommensurability of s, . This also mod- 
els a genèric case of random sequencing of successive s j . 
We have checked that different choice of Si, different or- 
dering of relative size Sj+i — Si does not alter the esscntial 
characteristics of the results. 

Despite the very simple analytic expression l|15|l of the 
scattering amplitude, as we increase N, the scattering 
quickly acquires "quantum chaotic" features [T^ |. which 
is characterized by Ericson fluctuation |15| . or the wild 






FIG. 3: Nearest neighbor spacing distribution P(s) obtained 
from the scaled energy levels of one-dimensional line limited to 
x — [-L, L] with scale invariant point interactions. Location 
parameters Sj are same, apart from the re-scaling, as in FIG. 
2. Dirichlet boundary condition is imposed at both ends, x = 
—L and L. The sòlid line represent the Wigner distribution 
P(s) = 7r/2 ■ s exp(— 7r/4 ■ s 2 ), and the dashed line the Poisson 
distribution P(s) — exp(— s). 



oscillation in scattering amplitudes caused by the over- 
lapping resonances. Because of the scale invariance of 
the each point interactions, the fluctuation appears in 
arbitrary energy scale. Clearly, this fluctuation is the 
result of interferences among multi-periodic oscillations 
with incommensurate freqüències, whose number of pe- 
riods 2 N ~ 1 proliferates very fast with increasing N. 

We next examine the statistical properties of energy 
level sequence calculated from our system with periòdic 
boundary condition. In FIG. 3, level spacing distribution 
P(s) for the system with several a are shown for N = 3, 
N = 5, and N — 7 cases from top to bottom as FIG3(a), 
FIG3(b) and FIG3(c). The distances Si are chosen to be 
= ^fpl : y/pi+i. The total length 



n+i 



L is set to be si : 2L = 1 



Eí=i Pi 



We have chosen 
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a = 27 to approximate disconnected large coupling limit, 
and a = 2 as the strong coupling limit, while, as an 
intermediate coupling, we chose the value a = 5. Thcsc 
graphs clearly show the approach of P(s) to the Wigner 
distribution (also known as GOE distribution), which is 
regarded as the quantum signature of classically chaotic 
dynamics |l6j | , at the strong coupling value as we increase 
the number of points N. The convergence appears to be 
fast as the Wigner-like level statistics already takes shape 
even with N = 3. We have confirmed, with numerical 
calculations up to N = 19, that very good convergence 
to the Wigner statistics is obtained with large N if the 
coupling paramcter a is properly rescaled to ensure the 
partial reflection/transmission of the system. 

We now discuss the implication of our findings in a 
broader context. The central result of this article is 
the generation of the "random" or irregular properties in 
quantum dynamics out of fully analytic quantum spectral 
functions obtained from a one-dimensional system. This 
type of properties are usually associated to the noninte- 
grable system. Since the classical counterpart of conser- 
vative one-dimensional system necessarily is integrable, 
the well-established correspondence between the chaotic 
classical dynamics and the random quantum dynamics 
seems to fail for our model. 

The clue to understand this seeming contradiction 
might be found in the singular nature of the high-energy 
limit of our system. Because of the special property of 
scale invariance present in Fülòp-Tsutsui point interac- 
tion, high energy limit, k —* oo does not bring the system 
to classical limit. As is well- known, high energy limit of 5 
and 5' interactions corresponds to the free pass and per- 
fect bouncing wall respectively, thus supplying two lcgit- 
imate deterministic classical limits of a point interaction. 



With scale invariant point interactions, we are presented 
with semi-transparent wall with finite penetration prob- 
ability at all energies. Therefore, if we were to identify 
the high energy limit as a classical limit, we are forecd 
to consider stochastic dynamics whose randomness orig- 
inates directly from the probabilistic nature of quantum 
mechanics itself. 

Irrespective to the problem of classical limit and corre- 
spondence, our analytical expressions, made possible by 
the scale invariance of Fülòp-Tsutsui point interaction, 
shed light on how irregular quantum dynamics emerge 
as the infinite-period limit of multi-periodic scattering 
matrices, just as chaotic classical dynamics emerge as 
the infinitc-period limit of multi-periodic motion. In this 
connection, it should be useful to consider a complc- 
mentary approach of trace-formula based analysis to our 
model. With appropriate modifications, existing semi- 
classical treatments of quantum graphs |Ï7^ appear ca- 
pable of handling the problem, and the comparison to 
the current approach should yicld further insight into the 
singular and irregular dynamics in quantum mechanics. 

Lastly, we emphasize the utility of scale invariant 
interaction in the studies of aspeets of one-dimensional 
systems other than the current example of stochastic 
dynamics, such as the spectral properties of regular 
lattice, in which analytical expressions (|16[) - H2U|) are 
expected to become very convenient. 
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